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The impacts of granular jets for both frictional and frictionless grains in two dimen- 
sions are numericaUy investigated. A dense flow with a dead zone emerges during the 
impact. From the two-dimensional simulation, we evaluate the equations of state and 
the constitutive equations of the flow. The asymptotic divergences of pressure and shear 
stress similar to the situation near the jamming transition appear for the frictionless 
case, while their exponents are smaller than those of the sheared granular systems, and 
are close to the extrapolation from the kinetic theoretical regime. In a similar manner 
to the jamming for frictional grains, the critical density decreases as the friction con- 
stant of grains increases. For bi-disperse systems, the effective friction constant defined 
as the ratio between shear stress and normal stress, monotonically increases from near 
zero, as the strain rate increases. On the other hand, the effective friction constant 
has two metastablc branches for mono-disperse systems because of the coexistence of a 
crystallized state and a liquid state. 

Subject Index J44, JOl, A56 



1. Introduction 

Non-equilibrium phenomena induced by impacts have been extensively studied in various 
contexts, such as nuclear reactions [1-3], nanotechnology [4, 5], water-bells [6, 7] and granular 
flows [8-16]. Crater morphology is studied via an impact process of a free-falling water drop 
or a grain onto a granular layer [9, 10], while a sinking grain produces a sand jet [11]. The 
impact of a granular jet on a target produces a sheet-like scattered pattern or a cone-like 
pattern, depending on the ratio of the target diameter and the jet diameter [8], which is also 
found in water-bell experiments with low surface tensions [6, 7]. 

Cheng et al. suggested that the fluid state of a granular jet after an impact is similar to 
the Quark Gluon Plasma(QGP), which behaves like a perfect fluid through their experiment 
[8] . Recently, we reported that the shear viscosity during the impact is well described by the 
kinetic theory of the granular gas [17-21], although the small shear stress observed in the 
experiment is reproduced through our three-dimensional (3D) simulation [12, 13]. Because 
the shear viscosity is not anomalous, the correspondence between a granular flow and QGP 
would be superficial. 

To discuss the fluid state of granular jets, we need to know the details of rheology of 
moderate dense granular flows. A typical situation of the study for a dense granular flow is 
the flow on an inclined plane [22-24]. Bagnold proposed the constitutive equation for dense 
granular flows that the shear stress is proportional to the square of the shear rate [22] , so 
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called Bagnold scaling, which is verified experimentally [23] and numerically [24, 25] under 
several conditions such as the flow down an inclined plane. Dense granular flows, however, 
have more variety of rheological constitutive equations for flows on inclined planes [26-30]. 
Conventional one would be the constitutive equation presented by Jop and coworkers [27], 
where the effective friction constant, defined as the ratio of shear stress to pressure, saturates 
from a static value at zero shear rate to a maximum value as the shear rate increases. The 
power-law friction law, which is also different from Bagnold's scaling for dense granular flows, 
is proposed via extensive simulations [31-34]. 

A granular system has rigidity above a critical value of density (j)j and does not have any 
rigidity below (j)j. This sudden change in the rigidity is known as the jamming transition 
[35-48]. The jamming is not only investigated in the system of grains, but also that of 
colloidal suspensions [49] or foams [50] . Here, (j)j decreases as the friction constant fip of 
grains increases. Moreover, it seems that there are two fictitious jamming points except for 
the true jamming point for finite [36]. Critical exponents of the divergence of the pressure 
and the shear viscosity near the transition are extensively discussed [36-48]. 

In this paper, we perform two dimensional (2D) simulations for the granular jet in terms 
of the discrete element method (DEM) [51]. Because grains are easily packed through the 
impact in 2D, the system would be near the jammed state. This paper is complementary to 
the previous 2D DEM study [16], and hard core simulations supplemented by the simulation 
of a perfect fluid model [15]. Indeed, although Huang et al. reported that the relevant role 
of the contact stress in a 2D granular jet, they were not interested in the critical behavior of 
jammed grains induced by the jet. Guttenberg suggested that the friction constant does not 
play a significant role, at least, in the scattering angle [14], while the effects on the jammed 
state induced by jets are unknown. 

This paper is organized as follows: After the introduction of our numerical model in Sec. 2, 
we analyze the profile of the local stress tensor, the area fraction and the granular tempera- 
ture, and discuss the rheology of the granular jets for the frictionless case in 2D to compare 
their behavior with the jamming transition for a bi-disperse frictionless case in Sec. 3. The 
effect of the friction constant is discussed in Sec. 4. In Sec. 5, our numerical results for a 
mono-disperse case are shown and the paper is concluded in Sec. 6. In the Appendix A, we 
comment on the artificial burst-like fiow in 2D, which appears in the case of large jip for soft 
grains. 

2. Model 

We adopt DEM to simulate the jet [51]. We mainly focus on bi-dispersed soft core particles of 
the diameter d and 0.8(i with the same mass m to avoid the crystallization. When the particle 
i at the position rj and the particle j at Yj are in contact, the normal force F^j is given by 
F^- = F^f + F^P""^ with Fj-f = kn{d - nj) and f/j"""^ = -y]n{^ij ■ rij), where = |rj - rj\ 
and gij = Vj — Vj with the velocity Vj of the particle i. The tangential force is given by 
Fj^j = mm{\Ffj\, fj,pFJj}sgn{Flj), sgn(x) = 1 for x > and sgn(x) = —1 for otherwise, Ffj = 
h^ij — Vt^ij with the tangential overlap Sjj and the tangential component of relative velocity 
5jj between i th and j th particles. We examine the value of fip from fj.p = 0.2 to 1.0. 
Here, we adopt parameters kn = 4.98 x lO^muQ/d^, r]n = 2.88uo/d, with the incident velocity 
uo for the frictionless case and fip = 0.2. The value /Up = 0.2 is close to the experimental 
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Fig. 1 (i) A typical snapshot of the simulation for the frictionless case with </>o = 0.90. 
Blue particles and green particles denote grains with diameter 0.8d and d, respectively and 
red particles are wall-particles, (ii) All of the corresponding contact forces among grains are 
visualized as blue colored arrows. 

value for nylon spheres [52]. We use /c„ = 1.99 x lO^muQ/d'^ , rjn = 5.75uo/d for Hp = 0.4, and 
kn = 7.96 X lO^mtig/d^, ?7n = W-15uo/d for fip = 1.0. This set of parameters implies that 
the duration time is tc = O.lOd/uo for the frictionless case and fip = 0.2, tc = 0.05d/uQ for 
/Up = 0.4 and tc = O.Old/uQ for = 1.0, the restitution coefficient for a normal impact is 
unchanged e = 0.75 for the frictionless case and for all fip. The reason why we adopt these 
parameters is that the artificial burst-like flow emerges, if we adopt the same tc for large /ip, 
as is shown in the Appendix A because of many overlaps among grains. For the tangential 
parameters, we choose kt = 0.2kn,r]t = 0.5T]n- We adopt the second-order Adams-Bashforth 
method for the time integration of Newton's equation with the time interval At = 0.02tc- 

An initial configuration is generated as follows: We prepare a triangular lattice with dis- 
tance between grains l.ld and remove particles randomly to reach the desired density. We 
control the initial area fraction (/>o/(/>ini = </'o before the impact as 0.30 < (l)o < 0.90 with the 
initial area fraction before the removal (/>ini = 0.612, 0.780 for the bi-disperse and the mono- 
disperse case, respectively, and 8,000 particles are used. We average numerical data over 
the time 180.0 < tu^/d < 300.0 after the impact. The initial granular temperature, which 
represents the fluctuation of particle motion, is zero. The wall consists of particles in one 
layer with the same diameter d and the same mass m, which are connected to each other 
and with their own initial positions via the spring and the dashpot with the spring constant 
kw = lO.Omug/d^ and the dashpot constant rjw = B.Orjn, respectively. 

A typical snapshot of our simulation and that of the contact force network are shown in 
Fig. 1 (i) and (ii), respectively. Blue, green and red particles denote grains with diameter 
0.8d and d, and wall-particles, respectively in Fig. 1 (i) and all of the corresponding contact 
force network among grains are visualized as blue colored arrows in Fig. 1 (ii). It is easily 
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found that the dense flow with the weh-developed contact force network emerges during the 
impact. 

We evaluate physical quantities near the wall in two regions: < x < 5d and 5d < x < lOd, 
where we call (a) and (b) layers in the followings, respectively. We use Rjct/d = 15.0 and 
i?tar/-Rjet = 2.2 with the target radius i?tar and the jet radius -Rjet- We adopt the Cartesian 
coordinate, where y = is chosen to be the jet axis, and divide the calculation region into the 
y direction y = —5Ay, — 4Ay, • • • , 0, • • • , 5Ay, with Ay = i?tar/5. Then we estimate physical 
quantities in the corresponding mesh region with kAy < y < (k + l)Ay [k = —5, —4, • • • ,4). 
Numerical data are averaged over ten initial configurations with the same (pQ and error bars 
in figures denote their variance. 

We calculate the stress tensor as in Ref . [53] . The microscopic definition of the stress tensor 
at r is given by 

i i<j 

where i and j are indices of particles, /u, = x, y, the contact force between i th and j th 
particles Fj^^ and denotes the summation over the particles denoted by i located at r. 
A is the are of each mesh at r and nj^(r) = f ^ — v^{t) with the mean velocity v^{t) in the 
mesh at r. 



3. Rheology of Granular Jets for the frictionless case 

In this section, our numerical results of granular jets, mainly for 2D frictionless cases are 
presented. The results for frictional grains will be reported in Sec. 4. In Sec. 3.1, the existence 
of the dead zone and the profile of the area fraction are discussed. After showing profiles 
of the stress tensor in (a) or (b) layer in Sec. 3.2, we evaluate the equation of state and 
constitutive equation to compare our system with the critical behavior of the jamming in 
Sec. 3.3 and 3.4, respectively. 
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Fig. 2 The profile of Tg and (/>/(/)max for the frictionless case with (pQ = 0.90 is shown in (i) 
and (ii), respectively. There exists the dead zone in (a) layer, which is denoted by red empty 
squares. Blue filled squares denote Tg in (b) layer, where the dead zone does not exist. Tg in 
3D, which is denoted by purple asterisks, is higher than those in 2D. Grains are well packed 
in 2D, while they are not in 3D. 
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3.1. Existence of the dead zone and the profile of the are fraction 

Chicago group demonstrated the existence of the dead zone near the target, where the motion 
of the grains is frozen, [14, 15]. In our 3D study, however, we confirm that there is no definite 
dead zone, while the average velocity of grains becomes smah near the target [13]. On the 
other hand, in this paper, we verify the existence of the frozen layer (a), while the motion 
is not frozen in (b) layer for 2D granular jets. The comparisons of the granular temperature 
Tg defined by Tg = Xljgc "''^i /-^-^c with the number of grains Nc in the mesh c in 2D and 
3D for the frictionless case are shown in Fig. 2 (i). Although Tg slightly drops at the center 
in 3D, it is much larger than those in 2D and those for perimeter regions in 3D. Note that 
the profile of the packing fractions divided by </>max with (/iniax = 7r/(2V3) ~ 0.907 in 2D and 
</'max = 0.740 in 3D are shown in Fig. 2 (ii) for the frictionless case. It is apparent that the 
granular fluid after the impact in 2D becomes denser than that in 3D. 

3.2. Profile of the stress tensor 

The profiles of the stress tensor for (a) and (b) layer of frictionless grains are shown in Fig. 3 
(i) and (ii), respectively. We stress that there exists a large normal stress difference between 
a XX and Uyy in each layer as in 3D case [12]. 
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Fig. 3 The profiles of the stress tensor in (a) and (b) layer for = 0.90 are shown in (i) 
and (ii), respectively. There exist large normal stress differences between a^x and cjj^y, and 
the shear stress is much smaller than the normal stress in (b) layer, though it is not so small 
in (a) layer. 



Ellowitz et al. suggested that the profile of the velocity and the pressure for the granular 
jet is reproducible from the simulation of a perfect fiuid [15] but our result may not support 
their claim. Indeed, the shear stress looks small but finite. Moreover, the large normal stress 
difference exists in both layers, which does not exist in the perfect fluid. We should note 
that they have not discussed the stress tensor itself in details, though they reproduce some 
similar feature of their hard core simulation. In addition, Huang et al. indicated the relevant 
role of the contact stress in their DEM simulation, which may be an indirect objection to 
the perfect fluidity of the jet flow [16]. 



5/15 



3.3. Equation of state 

Let us discuss the equation of state for the 2D granular jet impact. We estimate the strain rate 
D^y = {dvy/dx + dvJdy)/2 as dvy{l:^x/2,y)/dx ~ {vy{'i^x/A,y) - 7;j,(Ax/4, y))/(Ax/2), 
dvy{3Ax/2,y)/dx ~ {vy{7Ax/A,y) - Vy{5Ax/4:,y))/{Ax/2) and dv^{x,y)/dy ~ {v^{x,y + 
Ay/2) —Vx{x,y — Ay/2))/Ay. Since physical quantities are evaluated near the wall, the 
mesh < X < Ax is divided into < x < Ax/2 and Ax/2 < x < Ax, and Ax < x < 2Ax 
is divided into Ax < x < 3Ax/2 and 3Ax/2 < x < 2Ax to calculate dvy{Ax/2,y)/dx and 
dvy{3Ax/2,y)/dx. -/?tar <y < ^tar is divided into -llAy/2 <y < -9Ay/2, -9Ay/2 < 
y < -7Ay/2, •••9Ay/2 <y < llA?//2 to calculate dvx{x,y)/dy. 

We follow the analysis in Ref. [31]. Here, we introduce two dimensionless numbers consist- 
ing of pressure: It = ^/T^JPd? and Ig = Dxyyjm/ P with pressure P = [a^x + fj/i/)/2- We 
plot numerical data on (p vs It plane and (j) vs Ig plane, in Fig. 4 (i) and (ii), respectively. 
As can be seen in Fig. 4 (i), <j) in (b) layer has a little larger value than (j) in (a) layer at the 
same It, while all (p against Is are collapsed on a universal curve (Fig. 4 (ii)) 

We can fit the data by the equations 

= </>T-aT4^"" (2) 
cj, = <p^-aJ^J^', (3) 

with constants 4>T,aT,aT,4>s,as and a^. Fitting parameters are determined simultaneously 
by using Levenberg-Marquardt algorithm [54]. The obtained equations of states are written 
as 



P£ _ a?- 



Tg (0r - 

P a° 



mDly 



(4) 
(5) 



where the comparison of Eqs. (2) and (3) with numerical data for the frictionless case are 
shown in Fig. 4 (i) and (ii), respectively. From Eqs. (2) and (3) which suggest the pressure 
diverging at (pT or i;^^ , the granular particles are well packed with the fraction sufficiently 
close to the jamming point. 

The obtained parameters from our simulation are (pT = 0.858 it 0.006, aT = 0.980 ± 
0.1, ar = 1-15 ± 0.1, 0.834 ±0.001, = 3.94 ±0.5 and = 1.36 ± 0.05, where the 
error originates from the fitting. Note that the conventional jamming point (pj ~ 0.8425 
at which the pressure diverges is located between cps and (pT and close to (ps [48]. The 
asymptotic divergences of pressure for the frictionless case are described as 



Tg mD'^y 



(6) 



In a conventional picture based on the extrapolation of the kinetic theory, the divergence of 
the pressure is expected to originate from the divergence of the radial distribution function 
i.e. 

Pd? ^ 



1 = M<P) « — ^ (7) 



as — )• (pc, with the radial distribution function g{(p) and the critical density (pc [17]. In 
our case, the data are not far from Pd? /Tg ~ P/niD^y ~ (cpc — (p)~~^ expected from the 
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Fig. 4 Numerical data for the bi-disperse case of frictionless grains are plotted on (j) vs 
It plane (i) and (p vs Is plane (ii). Red and blue points denote data of (a) and (b) layer for 
several ^o, respectively. The corresponding solid lines in figures are fitting equations (2) and 
(3). 



conventional view based on the extrapolation of the kinetic theory, where Tg ~ mcPD'^y is 
assumed. 

On the other hand, Hatano demonstrated an elegant scaling law in the vicinity of where 
the corresponding exponents are estimated as = 2.8 and ax = 1.7 from his data of the 
jamming transition [45]. Otsuki and Hayakawa showed the phenomenological explanation of 
the critical behavior near (pj and they predicted Og = 4.0 and ot = 2.0 [41, 44]. It should be 
stressed that the critical scaling of the jamming transition is analyzed in the D^y — ?• limit, 
and the critical exponents strongly depend on the choice of the jamming point. Because 
the strain rate cannot be controlled in our setup, the jamming point is not clearly defined. 
Moreover, there are no data above the jamming transition in which the residual stress exists. 
Thus, our obtained exponents are smaller than those of the jamming transition for sheared 
granular particles. Although ax is close to the results of the kinetic theory, the analysis of 
ar and (pr may not be reliable, because the data for (a) and (b) layer are separated, due to 
the difference of the profile of Tg. 

Prom Eqs. (4) and (5), Tg and D^y are expected to satisfy 

Tg a'^^{^T-<pr-' 
The validity of Eq. (8) is verified in Fig. 5, which can be independent check of the scaling laws 
(4) and (5). From Fig. 5, Eq. (8) well reproduces the data for (p < (pg. This is a remarkable 
result which is in contrast to the conventional view Tg oc mcPDly [25]. 

3.4- Constitutive equation 

3.4- 1- Effective friction constant. Let us discuss Is dependence of effective friction con- 
stant ^* = —axy/P to obtain the constitutive equation. Numerical data for the frictionless 
case are shown in Fig. 6. The behavior of is conventionally described as 

^i*{Is) = /is + (Mmax - /is)(l + h/Is), (9) 

where /i* starts from a static value of ^s at zero shear rate and converges to a limiting value of 
//max at high We obtain ^s = 0.0153 ± 0.009, /i^ax = 0.521 ± 0.06 and Iq = 0.0820 ± 0.02 
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Fig. 5 The numerical data for \J md?Dxy/Tg with Eq. (8) are compared for (p < (j). 



for the frictionless case by fitting. Thus, f^*{Is) can be fitted by the conventional relation 
(9), which is denoted by a solid line in Fig. 6. It should be stressed that fj,s is close to zero. 
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Fig. 6 Numerical data for the bi-dispersed frictionless grains are plotted on fi* vs Is plane. 
Red and blue points denote data for (a) and (b) layer for several (po, respectively. All points 
are fitted into the phenomenological equation in Eqs. (9) and (10), where we cannot judge 
which equation is better from the data. 




Some researchers proposed a different constitutive equation called the power-law friction 

fi*=fis + bl^, (10) 



which also well reproduces numerical data [29-34], where /3 ranges from 0.28 to 1.0, depend- 
ing on the dimension, microscopic parameters and the friction constant of grains. Numerical 
data can be fitted by Eq. (10) within error bars, where we obtain 5 = 1.18 ±0.1 and 
/3 = 0.592 zb 0.03, assuming /i^ = for the frictionless case. The fitting result of Eq. (10) 
is denoted by a dotted line in Fig. 6. As can be seen in Fig. 6, there is no significant dif- 
ference between Eqs. (9) and (10). We, of course, cannot discuss the superiority of one of 
frictional laws from our simulation. 
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Fig. 7 Equation (12), which denotes the divergence of the shear stress, is independently 
checked for the bi-disperse frictionless case. Red and blue points denote data for (a) and (b) 
layer for several 0o, respectively. Eq. (12) well reproduces numerical results. Because axy 
itself in (b) layer are small, the error bars in (b) layer are larger than those in (a) layer. 



3.4-2. The asymptotic divergence of the shear stress. Let us discuss the asymptotic 
divergence of the shear stress: 



with an exponent fSg. By using the divergence of the pressure (5) and the power-law friction 



fi* oc /f , we obtain the constitutive equation for cj^ 



(12) 



which is checked independently against the numerical data (Fig. 7). The exponent is esti- 
mated to be /3s = (1 — (3/2)as — 0.96. Because we use the power-law friction /? > 0, the 
divergence of shear stress may be slightly weaker than that of P. Here, Bagnold's scaling is 
still satisfied even in the vicinity of the "jammed" density, which is contrast to the actual 
jamming transition [35-48]. 

Let us compare our observed critical behavior of the shear stress with the case of the 
jamming transitions, in details as well as that of the extrapolation of the kinetic theory. 
When we adopt the extrapolation of the kinetic theory, we have —axy/rriD'^y = rj/mDxy ~ 
?7* ~ (l)'^g{4>) ~ (f)'^/{(l)c — (j)) as (j) ^ 4>c, i-e. = 1-0, where Tg ~ md'^D'^y is used and dimen- 
sionless shear viscosity rj* = r]/rjQ is introduced with r/o = ^JmT^j^TuP and shear viscosity 
?7 = —(Jxy/Dxy The extrapolation from the kinetic regime by Garcia- Rojo et al. predicts 
that axy diverges at density different from Pd'^ /Tg and Ps = 1-0 [46]. Therefore, our analysis 
based on the power-law friction Eq. (10) predicts the results similar to Garcia- Rojo et al [46]: 
/3s = 1.0 and 0s < On the other hand, the exponents for the divergence of —axy/mD'^y 
at the jamming transition for sheared granular materials are estimated to be /3s — 2.6 from 
data in Ref. [45] and /3s = 4.0 in Ref. [41]. Thus, our corresponding exponent /3s = 0.96 is 
much smaller than those of the jamming transition for sheared granular systems and rather 
close to the result of the kinetic theory. 
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4. Effect of the friction constant 

We examine how the fluid state depends on the friction constant from the simulation for 
fip = 0.2,0.4 and 1.0. It is noteworthy that the separation between (a) and (b) layer exists 
for larger fip, even on (p vs Is plane. The results for fip = 0.2 and fip = 1.0 are shown in Fig. 
8 (i) and (ii), respectively. Because there are two branches on (/) vs Ig plane, we adopt Eq. (3) 
to fit the data in (a) or (b) layer, separately. Figure 9 denotes the critical densities and the 
exponents for each fip, where 4>s in both (a) and (b) layer slightly decrease as fip increases, 
and Q<j in (a) layer increases as fj,p increases, while it decreases in (b) layer. We note that 
the decrease of our critical densities (ps both in (a) and (b) layer are gentler than (pL of the 
jamming transition for sheared granular systems [36]. 
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Fig. 8 Fitting results of Eq. (3) for fip = 0.2 (i) and fip = 1.0 (ii). As fip becomes larger, 
data for (a) and (b) layer deviates from each other. 




Fig. 9 The fip dependence of the critical density (ps and the exponent in (i) and (ii), 
respectively, (ps in both (a) and (b) layer slightly decrease as fip increases. The purple solid 
line denotes the corresponding critical density of jamming for frictional granular particles 
(pL- The exponents in (a) layer increases as fip increases, while they decrease in (b) layer. 
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In contrast, the friction law is little affected by the friction of grains. The results for the 
friction law Eqs. (9) and (10) are shown in Fig. 10 (i) for /ip = 0.2 and (ii) for fip = 1.0, 
where the numerical data can be fitted by both Eqs. (9) and (10). We stress that the /j,* 
monotonically increases from near zero, as Is increases, even for the large fXp. 




Fig. 10 Numerical data for fip = 0.2 (i) and /ip = 1.0 (ii) can be fitted into Eqs. (9) and 
(10) within error bars, where we cannot judge which equations are better. The friction law is 
little affected by the friction of grains. It should be noted that the /i* monotonically increases 
from near zero, as the increment of Ig, even for the large /ip. 



5. Result for the mono-disperse case 

Here, we discuss the impact of granular jets in 2D for the mono-disperse case. A typical 
snapshot zoomed near the target is shown in Fig. 11, where grains are crystalized near the 
wall. The black solid lines in Fig. 11 (i) are drawn by hand to clarify the grain boundary 
between the crystallized region and the disordered region, where the boundary becomes a 
slip line. We also visualize all of the corresponding contact force network in Fig. 11 (ii). 

The notable difference of the mono-disperse cases from the bi-disperse cases appears in 
the friction law. We plot /i*(/s) for the mono-disperse case of frictionless grains, fip = 0.2 
and Hp = 1.0 in Fig. 12 (i), (ii) and (iii), respectively. First of all, n*{Is) for (a) layer and 
(b) layer cannot be fitted by a single curve, unlike the bi-disperse case. Judging from the 
snapshot (Fig. 11), grains, at least, in (a) layer are partially crystallized. Therefore, it is 
reasonable that the response of the crystallized region is different from that in disordered 
regions in (b) layer. 

The behavior of /i*(/s) in (a) layer, which are observed in both frictional and frictionless 
cases, can be understood as follows. Because of the crystallization, a grain is trapped in 
a crystallized region. However, as Ig increases, the grain can escape from the crystallized 
region. Thus, /i* decreases as Is increases. 

The macroscopic friction //*(/<j) for frictionless grains are different from that for frictional 
grains in (b) layer. The most remarkable difference between the frictionless and the frictional 
cases is the existence of peak of fi* at a small Ig for the frictionless case, while there is no such 
a peak for frictional cases. Because a frictional grain can roll over grains, grains easily form 
a cluster. Therefore, the boundary between such clusters becomes a slip line. Thus, fJ.*{Is) 
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Fig. 11 Typical snapshot of the simulation for the frictionless case with (pQ = 0.90 near 
the target. Yellow particles denote mono-disperse grains in the granular jet and red particles 
are wall-particles (i). The black solid lines are drawn by hand to clarify the grain boundary 
between the crystallized region and the disordered region, where the boundary becomes a 
slip line. Crystallization into a triangular lattice can be seen near the region enclosed by 
the black lines. All of the corresponding contact forces between grains are visualized as blue 
colored arrows in (ii). 



would be constant as Ig becomes smaller. On the other hand, because a frictionless grain 
cannot either roll over them nor slip, it is trapped in the crystallized region even for the large 
Ig. Thus, for frictional and frictionless cases exhibit different behaviors in (b) layer. 

However, we should stress that there exist two metastable branches for both frictionless and 
frictional cases. 
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Fig. 12 Numerical data for the frictionless case (i), fip = 0.2 (ii) and //p = 1.0 (iii) are 
plotted on n* vs Is plane. Red and blue points denote data for (a) and (b) layer for several 
(pQ, respectively. Unlike bi-disperse cases, fJ.*{Is) for (a) and (b) cannot be fitted into a single 
curve. in (i)-(iii) show similar behavior in (a) layer. However, in (b) layer, because 

frictionless grains cannot roll over the crystallized region, /i*(/s) for the frictionless case 
shows similar dependence on Ig to that in (a) layer, while the corresponding /i*(/s) for 
frictional cases do not. 
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6. Discussion and Conclusion 

We have performed two-dimensional simulations for the impact of a granular jet and dis- 
cussed its rheology. We confirmed the existence of the dead zone, as is reported in Ref. [15] 
at least in (a) layer, unlike our previous three-dimensional cases [12, 13]. There exists large 
normal stress difference, which has not been reported previously. The shear stress is much 
smaller than the normal stress, at least in (b) layer. We need to solve the inconsistency in 
(a) layer with Ref. [15]. 

We have analyzed the rheology of frictionless grains after the jet impact. We found that 
the pressure and the shear stress diverge with exponents similar to the extrapolations from 
the kinetic regime, and their exponents are smaller than those of the jamming transition for 
sheared granular systems. We adopted the power-law friction for /x* to obtain the critical 
exponent for axy/mD'^y. The discrepancy between our case and the jamming transition 
for sheared granular systems would originates from (i) our system cannot reach the true 
jamming transition and (ii) the uncontrollability of Dxy in our setup. The jamming point 
for a sheared system is located between (ps < <pJ < <pT and is close to (j)s Our analysis 
based on the power-law friction Eq. (10) is consistent with that by Garcia-Rojo et al [46], 
where axy diverges at the density different from Pd?' /Tg [46]. 

The effects of the friction of grains Hp have been discussed. Although Guttenberg [14] 
suggested that does not play a significant role, at least, in the scattering angle via the 
approximate hard-sphere method [55], we found that the existence of the friction affects 
rheology of granular fluids after the impact. The separation between (a) and (b) layer appears 
for larger /x^, even on <j) vs Ig plane. The critical fraction (pg decreases as ^p increases, which is 
similar to the behavior of critical fraction of jammed frictional grains (j)L- The corresponding 
exponent increases (decreases) as the increment of ^p in (a) layer ((b) layer). 

The effective friction constant lJ-*{Is) for the mono-disperse case has two branches because 
of the coexistence of the crystallized state and a liquid state. On the other hand, fJ.*{Is) for 
the bi-disperse case can be described by known constitutive equations for dense granular 
flow[27-34]. 
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Appendix A. On artificial burst-like fiows for the large fip case 

In this appendix, we comment on the artificial burst-like flow in 2D, which appears in the 
case of large /ip with softer grains than those in the text. After the impact of a jet composed 
of softer grains with large /ip, the burst occurs when a grain slips, because large tangential 
force can be accumulated before the slip of a grain. In Fig. Al, we show the time evolutions 
of Tg at —Ay < y < in (a) layer for (i) the frictionless case and Hp = 0.2, and (ii) /ip = 1.0 
with several stiffness, where Tg for the frictionless case and fip = 0.2 reaches the small steady 
values, while Tg raise many times after the impact for Hp = 1.0 with large tc, due to the slip 
events. As tc becomes smaller, the burst-like flows are suppressed. Thus, we use harder grains 
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for large ^p. Though there are a few small raises of Tg for the frictionless case, they are out 
of our averaging time. 



- fip = 1.0 tc = O.ld/uo 
■lip = 1.0 <c = 0.075d/uo. 
-Hp = 1.0 tc = 0.05d/uo 
■ lip = 1.0 = 0.025d/Mo- 




300 



Fig. Al The time evolution of Tg at —Ay < y < in (a) layer for the frictionless and 
fip = 0.2 (i) and Hp = 1.0 for the several stiffness of grains (ii). As tc becomes smaller, i.e. as 
grains become stiffer, the burst-like flows are suppressed. 
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